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Abstract 

Construction, in the framework of a Nonequilibrium Statistical Ensemble Formalism, of a Meso- 
scopic Hydro-Thermodynamics, that is, covering phenomena involving motion displaying variations 
short in space and fast in time -unrestricted values of Knudsen numbers-, is presented. In that way, 
it is provided an approach enabling for the coupling and simultaneous treatment of the kinetics 
and hydrodynamic levels of descriptions. It is based on a complete thermo-statistical approach in 
terms of the densities of matter and energy and their fluxes of all orders covering systems arbitrar- 
ily driven away from equilibrium. The set of coupled nonlinear integro-differential hydrodynamic 
equations is derived. They are the evolution equations of the Grad-like moments of all orders, de- 
rived from a generalized kinetic equation built in the framework of the Nonequilibrium Statistical 
Ensemble Formalism. For illustration, the case of a system of particles embedded in a fluid acting 
as a thermal bath is fully described. The resulting enormous set of coupled evolution equations 
is of unmanageable proportions, thus requiring in practice to introduce an appropriate description 
using the smallest possible number of variables. We have obtained a hierarchy of Maxwell times, 
which can be considered a kind of Bogoliubov's characteristic times in hydrodynamic and which 
have a particular relevance in the criteria for stablishing a contraction of description. 



I. INTRODUCTION 



It has been noticed that one of the compli- 
cated problems of the nonequilibrium theory^ 



of transport processes in dense gases and liq- 
uids is the fact that their kinetics and hydro- 



dynamics are intimately coupled, and must fluctuations 

'ye treated simultaneously (e.g., see Refs. fl- Hence; more advanced ap p roach es are re- 
6j). On this we may say that microscopic de- quired to m these restrictions . Consider first 
scriptions of hydrodynamics, that is, associ- near homogeneity> which implies validity in 
ated to a derivation of kinetic equations from the Umit of long wavelengths ( or wave num- 
classical or quantum mechanics and contain- ber Q approaching zero ). To go beyond it 
ing kinetic (transport) coefficients written in ig necessary to introduce a proper depen- 
terms of correlation functions, is a long stand- dence Qn Q vaM) in principle; for inter . 
mg traditional problem. An important as- mediate and s hort wavelengths (intermedi- 
pect is the derivation of constitutive laws ate to large wave numbers). In phenomeno- 
which express thermodynamic fluxes (or cur- logical theories this corre sponds to go from 
rents, as those of matter and energy) in c i ass i ca ] irreversible thermodynamics to ex- 
terms of appropriate thermodynamic forces tended irreversible thermodynamics fl-Q. 
(typically gradients of densities as those of TWs ig what hag been called generaUzed hy _ 
matter and energy). In their most general drodynarmcs, a question extensively debated 
form these laws are nonlocal in space and for decades by the Statistical Mechanics com- 
non-mstantaneous m time. A first kinetic- muri ity. Several approaches have been used, 
hydrodynamic approach can be considered to and a descript i on can be consulted in Chap- 
be the so-called classical (or Onsagerian) hy- ter 6 of the classical book on the subject 
drodynamics- it gives foundations to, for ex- by Boon and Yip [ 2 ]. Introduction of non- 
ample, the classical Fourier's and Fick's dif- local effects for describing mot i ns with in- 
fusion laws. But it works under quite restric- fluence of ever decreasing wave lengths, go- 
tive conditions, namely, local equilibrium; ing towards the very short limit, has been 
linear relations between fluxes and thermo- done in terms of expansions in increasing 
dynamic forces (meaning weak amplitudes in powerg of the wavenumber; which consists 
the motion) with Onsager's symmetry laws in what ig sometim es referred to as higher- 
holding; near homogeneous and static move- Qrder hydrodynamics (HOH). Attempts to 
ment (meaning that the motion can be well perform such expansions are the so-called 
described with basically Fourier components Burnett and super-Burnett approaches in the 
with long wavelengths and low frequencies, cage of mags motioU} and Guyer-Krumhansl 
and then involves only smooth variation in approach in the case of propagation of en- 
space and time); weak and rapidly regressing ergy (gee for example Refs . [ U ] and [12 ]). 



3- r 



An usual approach has been based on the ids under flow 17|-[l9|. Indeed, the nonlocal 
moments solution procedure of Boltzmann terms become specially important in minia- 
equation, as in the work of Hess jl3], us- turized devices at submicronic lengths {20 1, or 
ing a higher-order Chapman- Enskog solution in the design of stratospheric planes, which 
method. The Chapman- Enskog method pro- fly in rarefied gases in a density regime be- 
vides a solution to Boltzmann equation con- tween the independent particle description 
sisting of a series in powers of the Knudsen and the purely continuous description. An- 
number, K n , given by the ratio between the other particular problem to it related is the 
mean-free path of the particles and the scale one of obtaining the structures of shock waves 
of change (relevant wavelengths in the mo- in fluids for wide ranges of Mach numbers 
tion) of the hydrodynamic fields. Retaining [21[. Moreover, Burnett approximation of hy- 
the term linear in K n there follows Navier- drodynamics has been shown to provide sub- 
Stokes equation, the term in K\ introduces stantial improvement on many features of the 
Burnett-like contributions, and the higher- flow occurring in several problems in hydro- 
order ones [K\ and up) the super-Burnett dynamics, e.g. the case of Poiseuille flow [22 1 
contributions. and others [231 ] . 

A satisfactory development of a HOH be- The microscopic derivation of a HOH, to- 

ing also nonlinear and including fluctuations gether with the analysis of the validity of ex- 

1 a point in question. 



24| that for the case 



is highly desirable for covering a large class isting theories, is sti 

of hydrodynamic situations, and, besides It has been shown 

its own scientific interest, also for obtain- of Maxwellian molecules, whereas Navier- 

ing insights into present day technological- Stokes approximation yields equations which 

industrial processes. Also, we can mention are stable against small perturbations, this 

its fundamental relevance in Oceanography is not the case when are introduced Bur- 



and Meteorology ( e.g. jl^j ]. jl5|). and that nett contributions to the equations. It fol- 

it has been stated 16( that the idea of pro- lows that small perturbations to the solu- 

moting hydraulics by statistical inference is tions, which are periodic in the space vari- 

appealing because the complete information able with a wavelength smaller than a critical 

about phenomena in hydraulics seldom ex- length, are exponentially unstable. This fact 

ists; for example sediment transport and the has been called Bobylev's instability. More- 



more fundamental problem in fluid mechanics over, Karlin [25[ reconsidered the question 



of describing the velocity distribution in flu- looking for exact solutions to simplified mod- 



els: When a linearized ten-moments Grad engineering, soft-matter engineering, oil pro- 
method is used, and the Chapman- Enskog duction and petrochemistry, etc., which have 
method is applied to the model, does in an associated economic interest. MHT was 
fact there follow instabilities in the higher- initiated by the so-called Catalan School of 
order approximations. On the other hand, Thermodynamics as a large expansion of Ex- 
resorting to the Chapman- Enskog solution tended Irreversible Thermodynamics [301 ] . 
for linearized Grad ten-moment equations re- it can be noticed that nowadays two ap- 
summed exactly, solutions are obtained for pr0 aches appear to be the most favorable 
which the stability of higher-order hydrody- f or providing very satisfactory methods to 
namics, in various approximations, can be deal with hydrodynamics within an ample 
discussed. However, morerecently, Garcia- scope of nonequilibrium conditions. They are 



Colin and collaborators 



26| have extended Nonequilibrium Molecular Dynamics (NMD) 

Bobylev's analysis for the case of any inter- ^ anc } the kinetic theory based on the 

action potential, and have demonstrated that f ar reaching generalization of Gibbs' ensem- 

one can interpret the fact as to give a bound ble formalism, namely the Nonequilibrium 

for a Knudsen number above which the Bur- Statistical Ensemble Formalism (NESEF for 

nett equations are not valid, with no insta- short) [321]- [36(] . NMD is a computational 

bility involved. method created for modeling physical system 

Furthermore, inclusion of nonlinearity at the microscopic level, being a good tech- 

in the theory, in a Mesoscopic Hydro- nique to study the molecular behavior of sev- 

Thermodynamics (MHT for short and mean- eral physical processes. Together with the so- 

ing thermal physics of fluid continua), leads called Monte Carlo method are part of what 

to additional possible singularities, called hy- is known as numeric simulation methods |37|. 

drodynamic singularities, as, for example, We do here present an extensive derivation 

those described in Refs. j^-j^j. A satis- of a MHT on the basis of the kinetic theory 

factory construction of a MHT is highly de- founded on NESEF, quite appropriate to deal 

sirable for covering a large class of hydrody- with systems in far-from-equilibrium condi- 

namic situations obtaining an understanding tions involving the development of ultrafast 

of the physics involved from the microscopic relaxation processes, and displaying nonlin- 



level, and in the last instance gaining insights ear behavior leading, eventually, to instabili- 
into technological and industrial processes as ties and synergetic self-organization |38| 



40]. 



in, for instance, hydraulic engineering, food Within the framework of NESEF, but in a 



different approach to the one used here, an and in thermal equilibrium with an external 
alternative MHT was introduced by Zubarev reservoir. The general theory for the MHT 
and Tishchenko 41(, 42[. is built upon a generalization of Grad's mo- 



It may be noticed that the formalism can ments method for solution of, in this case, 

be extended to deal with the so-called non- a generalized kineticequation derived in the 

conventional hydrodynamics which is asso- context of NESEF [58]. 
dated to disordered media {43}, consisting 

in systems showing a complex structure of n THEORETICAL BACKGROUND 
a fractal-like (self-affine in average) charac- 
teristics, whose range of applicability and For building a nonlinear higher-order 



of physical interest is large 



44J. Fall on (generalized) hydro-thermodynamics on me- 



this problem the case of the distinctive be- chanical statistical basis, one needs to resort 
havior of polyatomic structures such as col- to a nonequilibrium statistical ensemble for- 
loidal particles, surfactant micelles, and poly- malism (NESEF) for open systems. Such for- 
mer and biopolymer (as DNA) in liquid solu- malism was developed step by step along the 
tions, which are classical examples of what is past century by a number of renowned scien- 
presently referred to as soft condensed matter tists whose contributions have been system- 



45|]. One particular case of apparently un- atized and generalized in a close structure, as 



usual behavior is the one associated to hydro- described in Refs. [32] to 
dynamic motion leading to a so-called non- According to theory, immediately after 

Fickian diffusion, described by a time evolu- the open system of N particles, in contact 

tion following a kind of fractional-power law with external sources and reservoirs, has been 



46|]. The nonequilibrium statistical thermo- driven out of equilibrium, the description 



mechanical aspects of complex systems in- Q f its state requires to introduce all its ob- 
cluding illustrations is reported elsewhere servables, their fluctuations and, eventually, 
a case invo lvin g hydro-thermodynamics higher-order variances. In most cases it suf- 



is given in Ref. 



48J- flees to take a reduced set of observables, 



In the present paper the conventional what implies in to have access to the so-called 

NESEF-based MHT is described in next Sec- one-particle (or single-particle), nj, and two- 

tion, accompanied with the study of a sys- particle, n 2 , dynamical operators for any sub- 

tem consisting of particles embedded in a set of the particles involved. This is so be- 

fluid which acts as a thermal bath at rest cause all observable quantities can be ex- 



5 



pressed at the microscopic mechanical level xhi(r, p) 

in terms of these operators (e.g. Refs. Q _ j Mp j ^ p ^ p/ . t) 

and [50|). ^ ^ 

On the basis of the construction of the x /i 2 (r, p, r', p') >. (4) 

nonequilibrium statistical operator 32[- 33], ' 



and taking into account the noted above fact 
that a complete description of the nonequi- 
librium state of the system follows from the 

. , , ni . . , l , . n , taken as normalized, as it should, with <b(t) 

knowledge of the single- and two-particle den- 



Hence, pit, 0) depends on the variables of the 
system of interest and p R on the variables 
of the thermal bath; both distributions are 



sity operators (or equivalently the associated 
reduced density matrices) which in classical 



ensuring the normalization of p, that is, 



mechanics are ^ = dT exp \ - d 3 rd 3 p F 1 (r, p; t) 



N J 

™i(r,p) = ^6(r-r i )5(p-p i ), W 

j=i xm(r,p 



n 2 (r,p;r,p') = ^ 8 (r - t s ) 8 (p - p,-) 

x 8 (r - r fe ) 8 (p 1 - p k ) , (2) x n 2 (r, p, r', p') \, (5) 



where r„ Pj are the coordinate and momen- and Fj and F% are the nonequ ilibrium ther- 

tum of the j-th particle, and r, p are called mod y nam i c variables conjugated to hi and 

field variables, the most complete nonequi- ^ mean [ U g that 

librium statistical distribution {4], 32]- (36] is 

the one built in terms of the auxiliary statis- ri -=,, s 

8 In Z(t) 

tical operator 8F 1 (r,p;t) = ~ Tr Wr ' P)/?(t ' ° )} ' (6) 

0) _^L_._ rr{B2(r , p;r , p , Wt , 0)} 

where p refers to the system of N particles (7) 

of mass m, and pR is the one associated to where lnZ(i) = (f){t) with Z(t) playing the 

a thermal bath of N R particles of mass M role of a nonequilibrium partition function 

taken in equilibrium at temperature T . The and 8 stands for functional derivative, in 

first one is given by complete analogy with the equilibrium case. 



p{t, 0) = exp j - <f>(t) - J d 3 rd 3 p F x (r, p; t) 



Moreover, dT is the element of volume in the 
phase space of the system, and for simplicity 



we have omitted to indicate the dependence 

on T of "hi, n 2 , p, 71, and that p R depends on H = H + H', (10) 

the phase point T R in the phase space of the 

where Ho is the kinetic energy operator and 

bath. 

We stress that p of Eq.(jl]) is not the sta- 
tistical operator of the nonequilibrium sys- = + * + (») 

tern, but an auxiliary one - called the in- con t a ins the internal interactions energy op- 

stantaneously frozen quasi-equilibrium" sta- erator fi u while accoun t s for the interac- 

tistical operator - that allows to built the tion of the system with the thermal bath, and 

proper nonequilibrium statistical operator, fy p is the energy operator ass0 ciated to the 

which needs to include historicity and irre- coup li n g of the system with external pumping 

versibility not present in p, hence it does not sources . Finally, e is an infinitesimal positive 

account for dissipative processes, and besides rea i num ber which is taken going to zero af- 

does not provide correct average values in the ter the traces in the ca l cu l a tion of averages 

calculation of transport coefficients and re- h ave Deen p er fo rmec i (it is present in a ker- 

sponse functions. ne j ^hat introduces irreversibility in the cal- 

We recall that the nonequiljhrium statis- culations, in a Krylov-Bogoliubov sense Q- 



tical operator is given by |4|, |32|-[36| y ). We stress that the second contribu- 

tion in the exponent in Eq. ([8]) accounts 

j^i eit'-t) 

for historicity and irreversible evolution from 



K E (t) =exp|lnp(t,0) - j 



^ | the initial time (taken in the remote past, 

df ^ n ^^ ^ [ X ^ t — > — oo, implying in adiabatic coupling of 



correlations, (see for example Ref. 32|), or 

with p(t, 0) of Eq. fl4]), and where 

alternatively, can be seen as the adiabatic 
p(t', t'-t)= exp U(t - t')c)p(t', 0), (9) coupling of the interactions responsible for 



relaxation processes [5l|]). Moreover we no- 
is the auxiliary operator carrying on the tice that the time derivative in Eq. ()8]) takes 
mechanical evolution of the system under care of the change in time of the thermody- 
Hamiltonian H (C is the Liouvillian operator namic state of the system (the first term in 
of the system meaning iCA = {A,H}). Usu- the argument i.e. t') and of the microscopic 
ally the system's Hamiltonian is separated mechanical evolution (second term in the ar- 
out into two terms, namely, gument, i.e. t' — t, see Eq. (Q), and that the 



initial value condition is 7Z £ (t ) = p(t ,0) for 
to — > — oo. 

The nonequilibrium thermodynamic space 
of states [52] associated to the basic dynamic 
variables hi and n 2 is composed by the one- 
particle and two-particle distribution func- 
tions 

fx ( r > P; = Tr {"1 ( r > P) Qe (t)} 

= Tr {n! (r,p)g(t,0)} (12) 

/a 0, P, r, p>; t) = Tr {f2 2 (r, p, r, p') g £ (t)} 

= Tr {n 2 (r,p,r',p') f?(t,0)} 

(13) 

where we indicate that for the basic vari- 
ables, and only for them, the average with 
the statistical operator g £ is equal to the one 
taken with the auxiliary operator 32j- 34]. 
The trace operation Tr is in this classical 
approach to be understood as an integration 
over phase space; hi and h 2 are functions on 
phase space and ~g and q £ functionals of these 
two. The knowledge of the two distribution 
functions f\ and / 2 allows to determine the 
value and evolution of any observable of the 
system as well as of response functions and 
transport coefficients. 

The knowledge of fi (r, p; t) implies com- 
plete information about the actual distribu- 
tion of particles, and therefore of the phys- 
ical properties of the system. Alternatively, 
knowing all the moments of the distribution 



allows to have a complete knowledge of its 
characteristics. A knowledge of some mo- 
ments is not sufficient to determine the distri- 
bution completely; it implies in only possess- 
ing partial knowledge of the characteristics 



of this distribution 



61| . Grad noticed that 



the question of the general solutions of the 
standard Boltzmann equation can be tackled 
along two distinct lines. One is to attempt 
to solve Boltzmann equation for the distri- 
bution fi itself in specific problems. Other 
is to obtain new phenomenological equations 
in an approach initiated b y M axwell 63] and 
continued by Grad [[64j, [66|] fit was called 
Grad's moments procedure 65]). These mo- 
ments produce quantities with a clear phys- 
ical meaning, namely, the densities of par- 
ticles and of energy and the fluxes of parti- 
cles of first and second order in a restricted 
fourteen-moments approach. 

In brief, the rth-order moment is the flux 
of order r 

lM(r,t) = J d 3 put\p)fi(v,p;t) (14) 

where ut^ is the r-rank tensor, s = n for par- 
ticle motion and s = h for energy motion, 



4 r] (p) 



—...Or — times)...— 

m m 



(15) 



4 r] (p) 



2m 



IP) 



that is, the tensorial product of r-times the 
vector p/m; for s = n, r = stands for the 



8 



density, r = 1 for the vector flux (or current), To proceed further, and give a clear illus- 

r = 2 for the flux of the flux which is related tration of the functioning of the theory, we 

to the pressure tensor field, and r > 2 for consider the case of a solution of N particles 

all the other higher-order fluxes. For s = h, of mass ra (the solute) in a fluid (the sol- 

r = stands for the density of energy and vent) of Nr particles of mass M. The former 

r > 1 for the respective fluxes. The density of is subjected to external forces - driving it out 

energy h(r,t) follows from the trace of rain , of equilibrium -, and the latter (the thermal 

namely bath) is taken in a steady state of constant 

p 2 £ f ^ equilibrium with an external reservoir at tem- 
perature T a . An analogous case, but at the 
The set composed by n(r,t), I n (r,t), q Uan t U m mechanical level, is the one of ear- 
lier, t) and h{r,t) is the one corresponding rierg embedded in the ionic lattice in doped 
to Grad's fourteen moments approach. Fi- or photoinjected semiconductors (see for ex- 



h{r,t) = J ^^/i(r, P ;t). (16) 



53| 



-Q). 



nally, the hydrodynamic equations are ample 

J^lM ( r , t) = j d 3 p U M (p) -|a (r, p; t) We write for the Hamiltonian 

(17) 

with r = 0, 1, 2, and where is to be in- H = H s + H R + W + H P , (18) 

troduced Eq.flZZD- Equations (ED consists of wher6) ^ firgt term Qn ^ ^ 
an enormous set of coupled nonlinear integro- 

N 2 , JV 
n 1 

differential equations. Evidently, it can be H s — -J- -\- _ y ([r^ — r fc |) (19) 

handled only in a contracted version, intro- ^ k 

ducing the hydrodynamics of order 0, 1, 2, is the Hamiltonian of the particles in the so- 

etc..., thus classified according to the last flux lute, consisting of their kinetic energy and 

that is retained in the contraction of the de- their pair interaction via a central force po- 

scription. Criteria for deciding the order of tential; the second term is 

the contraction must be established (see Ref. 



t 68l] ) . Hydrodynamics of order zero leads for n r p2 , n r 

H = I $ (IR — R I] 

the density to satisfy Fick' standard diffusion R 2-^ 2M 2 ' R M 

equation, the one of order one to Maxwell- (20) 

Cattaneo equation, and the other orders to which is the Hamiltonian of the particles in 

generalized Burnett and super-Burnett equa- the solvent, acting as a thermal bath, consist- 

tions. ing of their kinetic energy plus their pair in- 



teraction via a central force potential; more- iliary nonequilibrium statistical operator for 
over the particles embedded in the bath, is 

W = ± £ «-(fo-HJ) (21) ?(*.0)=-p{-^W-/^AFi(r,P;i) 

j=l ii=l 

is the interaction Hamiltonian of the par- 
ticles with the thermal bath, and Hp = 
^2iVext( T ii Pi> t) is t ne Hamiltonian associ- 
ated to the external force acting on the par- 
ticles of the system. where 

Under the stated condition that the bath 
is in constant thermal equilibrium with an ex- 
ternal reservoir at temperature T a , its macro- 
scopic state is characterized by the canonical 
distribution 



xwi (r,p) 



n^M)> ( 23 ) 
i=i 



Qj (t, 0) = exp [-fa (t) - j d 3 rd 3 p F 1 (r, p; t) 
x^r-r^p-p,)} (24) 



expf-Fi^p;*)}, (25) 



is a probability distribution for an individ- 

_ i f „ r - r I ual particle, with (i) and fa (t) ensuring the 

Q R = Z expt-/3 H R j (22) _ _ ' 

normalization conditions of ~q and ~g~j. 

where O = [k B T ]~ l and Z is the cor- The nonequilibrium equation of state 

responding partition function. The auxil- Q ^ ig ^ Qne rdating ^ variables 

iary nonequilibrium statistical operator of /t ( r> p; f ) Fl ( r , p; *) 
the whole system is the one of Eq. fl3]) and 

Eq. flU). But, for simplicity, considering a /i ( r ; P! *) = Tr {"i (r, p) g (t, 0)} 
dilute solution (large distance in average be- 
tween the particles) or that the potential V is 

screened (e.g., molecules in an ionized saline or 

solvent, e.g. [57| ), we can disregard the inffu- -Fi (r, p; £) = - In /i (r, p; £) . (26) 

ence of the two particle potential, and then Qn the other hand; the evolution equat i n for 

ignore n 2 , that is, taking F 2 = in Eq. @ re- f following from the NESEF-based kinetic 

r~i 

tammg only n x . In that case, we choose the theory; derived ag ghown in Ref ig the 

single-particle reduced density, n^p | T), generalized kinetic equation 
as the only relevant dynamical variable re- 

9 P( r >p;^) 

quired. Hence, g (t, 0), of Eq. (0}, the aux- Qj.h ( r ' P5 *) + ' V /i ( r ' Pi *) 

10 



+F (r, p; t) ■ Vp/i (r, p; t) — B (p) f x (r, p; t) of successive order, e.g. |62j). On this H. 

, T2i / s r „ „i „ / x Grad noticed that the question of the general 

-4 ] (p)©[VpV]/i(r,p;*) 

solutions of the standard Boltzmann kinetic 
- B 2 (p) [VpVp] fx (r, p;t) = 0, (27) e q Ua tion can be tackled along two distinct 

obtained in the Markovian approximation lines - 0ne is to attempt to solve Boltzmann 



34], 59], 60], where equation for the distribution f x itself in spe- 

cific problems. Other is to obtain new phe- 

P (r, p; t) _ p 

m m ' nomenological equations which generalize the 

usual (classical-Onsagerian) fluid dynamical 

F(r,p;t) = -V^(r,p;t)-B 1 (p) 

equations. The object is to show the transi- 
— Fjvl (r; t) , (29) ^ion f rom the Boltzmann equation in which 

with the explicit expressions, for the vecto- a state is § iven b ^ A ( r > P5 £ ) to the conven- 

rial quantities A a (p), Bx (p), F NL (r;t), the tional fluid description in which a state is 

second-rank tensors Af (p), B§ ] (p), and the S iven h V the densit y n ( r > the velocity field 

scalar B(p), together with a description of v ( r ^) > and the stress tensor T[2] (r,t), in 

the physical meaning of the several contribu- a su ffi cient generality to cover a broad class 

tions, are given in Ref. We also wrote of Problems. This approach was initiated by 



(it 



the symbol for full contraction of tensors. Maxwell |63| and continued by Grad M 

The distribution f x (r, p; t) that follows was called Grad's moments procedure) |65j. 
solving Eq. (l2T|) provides a complete informa- We do perform here an extensive gener- 
tion about the actual distribution of parti- alization of the moments procedure, consist- 
cles, and therefore of the physical properties ing into introducing the full set of moments 
of the system. Alternatively, if one knows of /i(r, p;t), of Eq. (127ft . in the variable p. 
all the moments of the distribution we do These moments produce quantities with a 
have a knowledge of its characteristics. A clear physical meaning, namely, the densi- 
knowledge of some moments is not sufficient ties of particles and of energy and their fluxes 
to determine the distribution completely; it of all order: the two vectorial fluxes, or cur- 
implies in only possessing partial knowledge rents, the tensorial fluxes, beginning with the 



on the characteristics of this distribution |61| second-order one which is the flux of the first- 

(this is related to Tchebychef 's procedure for order flux (the current of particles) which is 

obtaining characteristics of a probability dis- related to the pressure tensor, and all the 

tribution when we do possess the moments other higher-order fluxes. In that way we ob- 

11 



tain a quite generalized Mesoscopic Hydrody- where vfi is the /-rank tensor consisting of 
namics coupled to a N on- Equilibrium Ther- the tensorial product of I vectors u of Eq. 
modynamics, all together in the kinetic ap- (|32l) that is, 



proach provided by NESEF, as described in r . 

u w (p) 

next section. 



——...(I — times)...— . (35) 
m m m. 



We do have what can be called the family of 

III MESOSCOPIC HYDRO- hydrodynamical variables describing the ma- 

THERMODYNAMICS IN NESEF. terial motion, 1. e., the set 

u . „ {n(r,t), l n (r,t), {lj?(r,f)}} (36) 
Let us introduce, m the variable p, 

the moments of the distribution fi (r, p;i), with 1 — 2,3, .... 

namely On the other hand, we have the family of 



hydrodynamical variables describing the ther- 
mal motion, consisting of 
which is the density of particles; f 3 p 2 



n(r,t) = J d 3 pf\(r,p;t), (30) 



h(r,t) = J d 3 p^h(r,p;t), (37) 
I n (r, t) = I d 3 p u (p) ft (r, p; t) , (31) . o 



where 



/2 
d 3 P^~ ^/i(r,p;t), (38) 
2m m 

/p2 
d3p 2m~ Uil] (P) fl (r ' P; ^ ' 



with I n being the flux (current) of particles; 



(39) 

with I = 2, 3, ...; that is, in compact form 
jM(r, t) = J d 3 p u [2] (p) f x (r, p; t) (33) those in the set 

where = [uu] is the tensorial product of [h(r,t), I h (r,t), {4 /] (r,t)}| , (40) 

vectors u, with J„ being the second-order . 

which are, respectively, the density of en- 
flux (or flux of the first flux), a rank-2 tensor, 

ergy, its first vectorial flux (current), and 

which multiplied by the mass is related to the . . 

the higher-order tensorial fluxes. It can 

pressure tensor and r , 

be noticed that m this case of a parabolic 

jH(r,t) = y d 3 p u [l] (p) /i (r, p; t) (34) type energy-momentum dispersion relation, 

^(p) = p 2 1 2m, the set of Eq. (1401) is encom- 

are the higher-order fluxes of order I > 3 (the . . , , . T . 

passed m the previous one: In fact 

previous three of Eqs. ( )30|) . (l3Tj) and (|33|) 

are those for / = 0, 1 and 2 respectively), h(r,t) = — Tr{/^(r,t)} (41) 
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Ur,t) = jTr 2 {lW(r,t)} , (42) 

^/W(r,t) + V-^ +1 ](r,t) 
where Tr 2 stands for the contraction of two ut 



indexes, and, in general 



1 1 



lfhr,t) = jTr 2 {e%,t)}, (43) 

2 +4 il (r,t) + j| ] (r,t) + 4 1 i (r,t), (45) 



for all the other higher-order fluxes of energy. 
Hence, any flux of energy of order I is con- 
tained in the flux of matter of order I + 2. 
In what follows we concentrate the attention 
on the study of the hydrodynamic motion of 
the particles, with heat transport to be dealt 
with in a future communication in this series. 



where p(l, s) means that we must take a per- 
mutation of the first free index (1) with the 
s-th (s=l,2,3,...,/) free index of the Cartesian 



tensor 



when written 



F ext {r,t)lt 1] {r,t) 
in the indicial notation. Observe that the 
number of terms in the sum is given by the 
number of all permutations of I symbols in 
which I — 1 are repeated. All this ensures the 
IV. MHT EVOLUTION EQUATIONS correct symme try of this contribution, that 
IN NESEF a f u lly symmetrical tensor of order I. 

The several terms on the right of Eq. (1451) 



We proceed now to the derivation of the 

are, 

MHT equations, that is, the equations of evo- 

J= ext (r, t) = -VV ext (r; t) , (46) 

lutions for the basic macrovariables of the 

family of material motion, i. e. those in set that is the applied external force, created by 
of Eq. (|36|) . the action of the potential V ext , and the terms 

Let us consider the flux of order I (I = Ai (p), B x (p) and F NL (r;t) present in Eqs. 
0, 1, 2, ...); its evolution equation is (JSHD and ([29]) have respectively been incorpo- 

f) r f) rated into Jr\ jfj and jfj L , which are, 



dt 



ij?(r,f) = J d 3 p U N (p)^/i(r,p;t). 

(44) 



T m, \ n R y/2n Mp o x - (Q)\ 2 

Using Eq. f lzH) . but excluding a dependence J| J (r,t) = — > x 

1/ m~ t * 



on p of the external force, after lengthy but 



V m 2 ^ Q 



straightforward calculations we arrive to the f d 3 p^^ |^QQu'' _2 '(p) 

general set of coupled equations for the den- ' p 



sity, 1 = 0, the current, 1 = 1, and all the 
other higher-order fluxes, I > 2, given by 
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x exp 



Q 2 



h (r, p; *) 



n R (Mf3 f\^\^(Q)\ z ( 1 1 



V v^m 2 



Q 



d 3 ^ [Qu^p) 



M m 

Q p 



X 



2?i 



(53) 



Observe that the notation 



x exp 



/i(r, P ;t), (47) 



with a = M(3 /2m 2 , if> (Q) is the Fourier 

transform of the interaction potential be- = |^Q— ... (I — 1) times (54) 
tween the particles and the thermal bath, i. 
e., Eq.f l2"Tj) and according to Eq. f[3"5"j) 



means that one has to sum all permutations 
of the vectors in order to ensure that the 



QQ— ... (I — 2) times 

m m 



censor has the same symmetry of the tensor 
[i] 



QQu^p) 

(48) I* 1 on ^ e nanc i side of Eq. ( 145]) . 
r ,i 1 r r , ,, Next, making a Taylor series expansion of 

Jnl M) = — ^p(M) Fjvl ] (M)L 

(49) 



le exponential in both contributions in Eq. 



(T47]) and in Eq. i. e. 



F Wi (r; t) = / dV / d s p'G NL (r' - r, p' 



Gjvl (r' - r, p' 



/i(r',p';0: 



(50) 



exp 



-"'o- p 



£rilQPa] 0p |a] 



fc=0 



A:! 



V Q 



x{zF(Q,p') 



where QW stands for 

QW = [Q... (/- times) ...Q] 



(55) 



(56) 



+ 



2vr 



1/2 



7T Q p' 

m Q 
x exp [z'Q 



xexp 



r — r 



|and we recall that © stands for full contrac- 
tion of the two tensors of rank 2k. Then us- 
(51) 

ing Eq. (|53p for F(p, Q), we can rewrite Eqs. 
( 14"T|) . (I49p and (|52|) in a closed form in terms 
of all the fluxes, namely 



V m 2 ^ Q 2 
Q 



xQ-V / ^^[QuI'-^p)] F(p,Q) 



jW(r,t)=^p / [Air 2] ©^- 2l (r,t)] 



k=0 



x /i(r, p;t) 



(52) 



oo I 



££p(M) [Al 2fc+2l 0/f +/] (r,t 



fc=0 s=l 



where 



(57) 
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oo I 



4 q («-.*)=EEp( 1 ' a ) 



k=0 8=1 



A 



[2fe+2] 



+ 



0Vr w, (r,4 (5 

-. I OO „ 

4(M)=--Ep(M)E/ iV 

8=1 fc=0 

AgS 1 l(r'-r)0JP*V,*)]/r i] (r,*) 



(59) 



where 



^Ekm+e^ 1 ' 



s=2 



r=3 



second involves local couplings with the dif- 
ferent fluxes. On the other hand J^ L is non- 
linear (bilinear) in the fluxes accounting for 
nonlocal correlations involving all of them. 
Next, we reorganize these expressions setting 
into evidence the contributions that contain 
the neighboring fluxes to the one of order Z, 
namely In~^ (r, t) and In 1 ' (r, t), that is, aris- 
ing out of the terms with k = and k — 1 
in the sum over k in Eqs. ( 157|) and (158]) ; to 
obtain Eq. ( H5]) in the form, 

1 ' 

-Ep(l.a)[^«t(r,t)/ B p - 11 (r ) <)] 



+ E E KV;2,s) (60) 

r=3 s=r+l 

is an operator involving the set of permuta- 
tions that ensures the proper symmetry of the 
tensor on which it acts. Here the operation 
A^l © A^l indicates the contraction of some 
indexes in order to give the right tensorial 
order of the equation. For example in Eqs. 
([57]) and ([58]) it indicates the contraction of 
(r + s — l)/2 indexes as to produce a tensor 
of rank I. The several tensorial kinetic coef- 
ficients are given in Appendix A: Eqs. flAlj) 
to Eq. (jMp . 

We can see that the expressions for jfi 
and j| ] are linear in the hydrodynamic basic 
variables of the set of Eq. 03 6 p with tensorial 
coefficients A' r l The first one contains con- 
tributions of a relaxation character, and the 



8=1 



~ df 1 ^ (r, t) + a io J>(1, s) [V (r, *)] 

s=l 

+2la L1 V ■ I^(r,t)+b TO {p t [l^M)] } 

+jH L (T,t) + S®(v,t). (61) 

The last term on the right of Eq. ( 16 ip is 
given by, 



5f(r,fHi-i^, 
m 



+3la Tl -ll\r,t) 
m 



+ a «-Ert 1 '*) i v/ r i] (r,t) 

1 

+ ^ ] (r,t), 



8=1 



(62) 



where Rn(r,t), given in Appendix A, con- 
tains the contributions of the fluxes of order 
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higher than I + 2; and the kinetic coefficients the present communication (case of a dilute 

a r i, dLo, a>Li, Ko, Ki are given in Appendix B. solution). However, as already noticed, such 

The first three contributions on the right of set of equations is intractable, and, of course, 

Eq. (J52"|) are associated to the fluxes of energy we need to look in each case on how to find 

of orders I — 2, l — l and I, terms that can be the best description using the smallest pos- 

considered consisting of thermo-striction ef- sible numbers of variables. In other words 

fects and which, then, couple these equations to introduce an appropriate - for each case 

with the set of kinetic equations describing - contraction of description: this contraction 

the movement of energy, i. e. the hydro- implies in retaining the information consid- 

dynamical variables of the set of Eq. (flOl) . ered as relevant for the problem in hands, and 



However it can be noticed the already men- to disregard nonrelevant information 67] . 
tioned fact that the fluxes of energy can be Elsewhere 



we have discussed the ques- 



given in terms of those of particles, namely, tion of the contraction of description (re- 
if = (m/2)Tr 2 |/n +2] |, [ cf. Eqs. flUJ) to duction of the dimensions of the nonequilib- 



(1431)]. Moreover, rium thermodynamic space of states). As 

shown, a criterion for justifying the different 

Of 1 = l[\a TO \ + {l-l)\b Tl \], (63) 

levels of contraction is derived: It depends 
with 6i playing the role of a Maxwell- on the ran g e of wavelengths and frequen- 
characteristic time for the l-th flux. cies which are relevant for the characteriza- 
We stress that I = corresponds to the tion > in terms of normal modes, of the hydro- 
density n(r,t), I = 1 to the first flux (cur- thermodynamic motion in the nonequilib- 
rent) I n (r, t), I = 2 to the second flux I l n ] (r, t) rium °P en system. Maxwell times have a par- 
which multiplied by the mass m is related to ticular relevance then we proceed to analyze 
the pressure tensor P^(r, £)], / > 2 to the them, 
other higher-order fluxes. Hence, Eq. ( IBTj) 

represents the coupled set of evolution equa- y THE HIERARCHY OF MAXWELL 

tions involving the density and all its fluxes in CHARACTERISTIC TIMES 

its most general form. It must be noticed that 

it is linear in the basic variables; no approx- Let us now analyze Maxwell characteris- 



imation has been introduced. Nonlinearities tics times of Eq. (1631) to show that they 
should arise out of the interparticle interac- follow a hierarchy of values. First we note 



tion inclusion which we have disregarded in that, taking into account Eq. flB7j) into Eq. 
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( |63|) and the definition in Eq. (IBlj) we can the third flux leads to the results that 83 — 
write 8i/3 and 6*3 ~ 56 l i/21 for the Brownian and 



8- 1 = I 



0i (64) 



Lorentz particles respectively. For any / we 
do have approximately: 
what tell us that any characteristic time for 1) for the Brownian particle (m/M » 1) 



5m+M v ; 



81 ~ (68) 



/ > 2 is proportional to the one of Z = 1, 
that is the one for the first flux which multi- 
plied by the mass m is the linear momentum 2) for the Lorentz particle (m/M << 1) 



density and then all are proportional to the q ^ ^ ,^ _|_ n n $ (gg) 

linear momentum relaxation time. On the 

other hand we do have that or ^ ~ 56 ^ P for lar S e 1 

Moreover according to Eq. f )65|) as the or- 



I 5(l + x)+/-l 



der of flux largely increases its characteristic 



q 1 _|_ 1 B(l-\-x)+l ' Maxwell time approaches zero, and 8i + i/8 



1 ^ 



r, 100 j -i /„, j 1, with both being practically null. In Fig. 

for / = 1,2,6,... and where x = m/M, and o *■ j & 

4.1 41 1 • 1 it is displayed the ratio of characteristic 

then the ordering sequence J 

Maxwell times, for flux of order I with the 
momentum relaxation time, ClS db function of 
81 > 8 2 > 83 > ... > 81 > 0j+i > .... (66) m /M, i. e ., the quotient of the masses m of 

...... . . . , the system and M of the thermal bath. 

is verified wfiicfi can be considered to rep- 
resent a kind of Bogoliubov's hierarchy of 



characteristic times 



50] in generalized hydro- 



dynamics, and we can see that 81 — > as 
I — > 00. Moreover according to Eq. fl64|) it 
follows that 

_ 5(l + s) 

Bl ~ i[Hi + x) + (i-i)} dl - (67) 

Comparing with the second flux {1 = 2), Figure 1: The quotient between several 

the one related to the pressure tensor, it fol- Maxwell characteristic times and the one of 

lows that for the Brownian particles (x » 1) the first flux as a function of m/M. 

82 — 8\/2 and for the Lorentz particles Observe in figure-1 that the quotient of 

(x « 1) 82 — 5#i/12. A comparison with masses has little effect on 81/81 but it varies 

17 



significantly with the order i of the flux and the system. 



681 ]. it can be con- 



plays a particular relevance in the criteria for As shown elsewhere 

stablishing the contraction of description dis- jectured a general contraction criterion, 

cussed below. namely, a contraction of order r (meaning 

keeping the densities and their fluxes up 

to order r) can be introduced, once we can 
VI. ON THE CRITERION OF CON- ' 

show that in the spectrum of wavelengths, 

TRACTED DESCRIPTION AND AN 

which characterize the motion, predominate 

APPLICATION 

those larger than a "frontier" one, A^ rr+1 ^ = 
Returning to the question of the con- v 2 6 r 6 r+1 , where v is of the order of the ther- 



tracted description it can be shown [fig] that mal velocity and 9 r and 9 r+1 the correspond- 

a truncation criterion can be derived, which in S Maxwell times. We shall try to illustrate 

rests on the characteristics of the hydrody- the matter using a contraction of order 2. 

namic motion that develops under the given Let us ^ TS ^ write down the equations of 

experimental procedure evolution, whose general expression is given 

Since inclusion of higher and higher order in Ec l- TO corresponding to the density and 

fluxes implies in describing a motion involv- its ^ nxes °f a11 order, for £ = 0: the density, 

ing increasing Knudsen numbers per hydro- for ^ = 1: the first flux of the density, ^ = 2: 

dynamic mode (that is governed by smaller the flux of the first flux which multiplied by 

and smaller wavelengths - larger and larger m * s tne P ressur e tensor field, and i = 3: the 

wavenumbers - accompanied by higher and A ux °f the pressure. We do have, respectively, 

higher frequencies), in a qualitative manner, — n(r t) + V I (r t) — (70) 
we can say that, general "thumb rule," 
the criterion indicates that a more and more 

restricted contraction can be used when larger _^_ T / ,\ . ^ T [2], ^_ n ( r ' ^) -r / +\ 
and larger are the prevalent wavelengths in 

the motion. Therefore, in simpler words, ^-n(Y,t) + aL \/n(r,t) 

when the motion becomes more and more + 2oliV • i^ 2 '(r, t) + S n (r, t) , (71) 

smooth in space and time, the more reduced 

can be the dimension of the basic macrovari- d [2l . r 3 i 

—I^{r,t) + V-r n i {r,t) 

ables space to be used for the description of 

the nonequilibrium thermodynamic state of — ~ {[^ ( r ^)^-n(^,t)] + [ I n (r, t) J-(r,t)]} 
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-9^I^(r,t)+a Lo {VI n (r,t) + [VI n (r, t)f} 

+ 4 a L1 V-I [ ^(r,t) + b TO n(r,t)l^+S^ (r,t) , 

(72) 

|^r,*) + V-/M(r,t)= 

s=l 

3 

+a LO $>(M) [VJ?(M)] 

+6 a L iV-4 4] (r,i) 

+ b TO {p 3 [l®I n (r,t)]} + S®(r,t), (73) 

where J- (r, t), given in Eq. ( IClj) . and the ex- 
pressions for S n (r, t), 5n (r, i) and Sl^ (r, t) 
are given in Appendix C. The Maxwell times 
Qy , O^ 1 and 9$ are obtained from Eq. (|63|) 
respectively for £ = 1,2 and 3. Moreover, as 
noticed, if we multiply Eq. ( J72l) by the mass 
m, we do have an equation for the pressure 
field tensor 

pW(r,t) = mI [ *\r,t), (74) 

composed of the hydrostatic contribution 
(the diagonal terms) and the shear stress (the 
non-diagonal terms) and the convective pres- 
sure (cf. Eq. (l9Tj) presented later on, but 
where the shear contributions have been dis- 
carded). We also mention that taking into 
account Eq. (1881) below relating I n (r, t) with 
the barycentric velocity v(r, t), Eq. ( I7T|) 
can be transformed in an evolution equation 
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for the latter to obtain a generalized Navier- 
Stokes equation (future publication). 

Let us now, for illustration, consider the 
case when we can perform a truncation in a 
third order, that is, to consider as basic vari- 
ables n(r, t) its flux I n (r, t) and and the pres- 
sure tensor m/jf (r,t). In this contracted de- 
scription we consider Eqs. (!70|) . ( ITTj) and ( I72|) 
but with the further restrictions in Eqs. ( I7T1) 
and (1721) of neglecting: 1. the shear stress 
contribution, more precisely introducing the 
trace of the pressure tensor which is propor- 
tional to the energy density h(r,t), that is, 

Tr (P [2] (r,t)} = 2h(r,t) (75) 

where convective pressure has been discon- 
sidered, cf. Eq. (159]) . 2. the terms with coef- 
ficients cll whose origin is in self-energy cor- 
rections, which simply would renormalize the 
kinetic coefficients, and 3. the terms S„ (r, t) 
and Sn (r, t) which contain the energy den- 
sity h(r,t) and its flux l h (r,t) thus, disre- 
garding thermo-striction effects. 

The evolution equations for the chosen 
hydrodynamic variables, n(r,t), l n (r,t) and 
7n (r, t), in the conditions above stated take 
the form 

^n(r,t) + V-I„(r,t) = 0, (76) 
|l n (r,t) + V-/L 2 '(r,t)=^^(r,t) 



p®(r,t) + V ■ I®(r,t)=b TO n(r,t)lM 



(77) Eq.( !74|) ]. For that purpose we consider con- 
ditions such that the pressure is changing 
slowly in time (8 2 dln ] /dt « 1$ or u 2 « 
1 along the motion), and then from Eq.( 178|) 
follows that 



m 



{[^(r,t)I n (r,t)]+ [I n (r,t) F(r,t))}. 

(78; 



V-i?(r,t) =b ro n(r,t)l 
-« ] (r,t) 



[2] 



Deriving in time Eq. (J76D and, next, in + L {[p ( T t )I n (r,t)]+ [I n (r, t)F(r, t)}} 

. m 



the result inserting <9I n (r, t)/dt, using Eq. 



(81) 



(J77D we arrive at the hyperbolic (Maxwell- and uging thig regult in Eq it foUows 

Cattaneo type) evolution equation for n(r, t) ^ e e q Ua ti n 



S n(r ' t) + ^ n(r ' t)=V ' V ' / " 1(r ' t) 



d 3 , . 
^ (r ' t) + 



1 1 

9i 2 



d 2 , x 
^(r,t) 



n(r, t) 



^M) 



(79) 



1 9 



0i0 2 dt 



r,t) = iv-V-jP(r,t) 



Furthermore, deriving in time this Eq. ([79 
and using Eq. flTHj) it follows that 

a 3 



-v- 



n(r, t) 



^(r,t) 



dt 



n(r, i) 



dt 3 



n(r,t) + 



1 1 

01 02 



<9f 



^M) 

(82) - 



2 ™M) 



+ ^^(r,t) = &™V 2 n(r,t) 



To close this Eq. (J82J) it is necessary to 
evaluate In which, we recall, is given by 



-V-V-V-4 3 l(r,t)--V- 

t>2 



1 



n(r, £) 



^M) 



J^(r,t)= / d 3 p Pj/^r.pjt). (83) 
mm 



n(r,t) 



m 



We resort now to the use of Eq. (1251) and for 
namely 



— V-V-{[^(r, t) I n (r, 0]+ [I„(r, t)jF (r, £)]} i*i(r, p, t) we use an expansion in variable p, 
m 



The divergence of the third-order flux, -^i( r 'P>^) Fi n (r,t) + 
V ■ /jf , in terms of the basic variables is 
evaluated on the basis of Eq. (l75|) . the evolu- 
tion equation for the second-order flux, which 
we recall, related to the pressure tensor [cf. 



0Fi(r,p,t) 



dp 



1 d 2 F 1 {r,p,t)\ tf_ 
2m 



P 



dp 2 



| 1 dPFtfafrt) 
2 dpdp 













PP 


+ 










(84) 
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where lower index nought indicates that the (called quasitemperature 69J, 70j), as well 



PP 



is the as 



derivative is taken at p = 0, and 
traceless part of the tensor. We rewrite iq in 
the form 



I h (r,f) = -n(r,t)0 _1 (r,t)v(r,t) 



in 



iq(r, p, t) ~ cp n (r, t) + F n (r, t) ■ p/m 



f^(r,t)v 2 (r,t)v(r,t), (90) 



and 



that is, keeping terms up to second order in 
p. This is consistent with the contracted de- 



I^(r,t)=-r\r,t)n(r,t)l 
m 



[2] 



+ n(r,t)v(r,t)v(r,t). (91) 



scription we used, and of disregarding the introducing Eq. flUD in Eq. flgS]), if 



n{r,t) 
m 



shear stress, and where </? n , F n and F h are y. 
the nonequilibrium thermodynamic variables we nna iiy arrive at 
conjugated to the density, the flux and the 

[21 

trace of I n which is proportional to the en- 
ergy. Moreover, introducing the alternative 
forms 

where 



» fc|v- 



n(r,t) 
m 



9 1 9^n{v,t) + [9 1 + 9 2 ] ^n(v,t) 

+ ^n(r,t) = -V-j n (r,t), (92) 



F n (r,t) = -m/3(r,t)v(r,t), (86) = -T>W(r,t)-Vn (r,t)-n (r,t)V (r,t) 

(93) 

plays the role of a generalized flux with at 
the right being present a generalized thermo- 
dynamic force, and where 

^l' 2 l + [vM)v(r,t)] 

m 

(94) ■ 

is playing the role of a generalized diffusion 
tensor (composed of two parts, a first one of 
^ ( r ^) — ~o^ r i^n\ r ^)} ~ y n ( r 5^0 P ( r 't) thermal origin and a second associated to the 



F h (r,t)=P(r,t), (87) 
it follows that 

I n (r,t) =n(r,t)v(r,t), (88) 

defining the barycentric velocity v(r, t) . ( r >0 — ^i 
From Eqs. (|30|), (151]) and (155]) there follows 
that 



m 



+ yn(r,t)v 2 (r,t), 



i9) 



drift of the material) and 



where we can write [3 1 (r,t) = ksT*(r,t) 
introducing a nonequilibrium temperature 



V(r,t) = 0! 



k B T*(r,t) 



ni 
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01 



V .[vv]-i^M) 

m 



(95) 



composed of three terms, one of thermal ori- w (r) — — exp {— r 2 /r 2 } , (98) 



gin another coming from the drifting move- 
ment and a third from the applied force. 

Moreover, in the steady state (dn/dt = 
and then j = 0) there follows that the density 
satisfies the equation 



with the open parameters U and r Q ( r Q is a 
length scale playing the role of a range length 
and U/r being the interaction strength ) . 
It has been noticed that this kind of po- 
tential belongs to the class of interactions 
V [2 \y) ■ Vn (r) = -n (r) V (r) . (96) which do not diverge at the origin, i.e., are 

bounded. They are potentials corresponding 
to effective interactions between the centers 
of mass of soft, flexible macromolecules such 



Returning to Eq. (|92|) . its Fourier trans- 
form reads 



n 

[-u 2 9 l e 2 - iu (6>i + 2 ) + 1] n(Q, u) as polymer chains [72|, dendrimers |73j, and 



Q. (ni \ / n -\ others. The centers of mass of two macro- 

•Jn(Q,w), (97) 

molecules can coincide without violation of 
which give us an illustration on the criterion tne excluded volume conditions .hence im- 



of contraction of description: 1. In conditions plying in a bounded interaction 74]. Several 
such that ujH^ « 1, the term with third studies of this potential can be consulted, for 
time derivative can be neglected and the evo- example, in Refs. 75] 



lution equation acquires the form of a gener- It can be not i ce d that GCM of Eq.® 

alized hyperbolic Maxwell-Cattaneo one. 2. roug hl y mimics a hard sphere potential with 

If further u {9 l + 2 ) « 1 the second time radius r ^ and that in the limit of To going 

derivative also can be neglected and we are to zero goes over a contac t potential U 5 (r) . 
left with a generalized parabolic diffusion- like p our j er transform is 

equation, and in that way there follows a f 1 1 

chain of increasing contractions of description ^ (Q) — ^p2^ T ° 6X ^ \ _ 4 T °^ J ' ^ 
of the hydrodynamic motion. 

In terms of these results we find that 

Finally, to perform numerical calculations 

iiii-i 1/1-1 V^" n R M ( M\ „o /2rr2 

and analyze the results we introduce a central 6, = — = — — 1 H p„' U , 

Q^2M l l 2 m\ m J ° 

force interaction between particles in the sys- (100) 

tern with those in the bath of the Gaussian what tells us that the momentum relaxation 

form, called the Gaussian core model (GCM) time becomes very large for the Brownian 
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particle when m » M, and very small for of macrovariables is given we can obtain the 

the Lorentz particle when m « M. Fur- nonlinear hydrodynamic equations, which are 

thermore, that B\ increases with the power the average value over the nonequilibrium 

3/2 of the temperature T Q and, as expected, ensemble of Hamilton equations of motion 

with the reciprocal of the density of scatter- (in the classical level or Heisenberg equa- 

ing centers. tions at the quantum level) of the basic mi- 

Moreover, it can also be noticed that in crovariables (mechanical observables) . Once 

the limit of a contact potential (r G — > 0), the all these hydrodynamical variables, cf. Eqs. 

quantity (1501) to (131]) and (f3"7|) to (1551). involve the sin- 
gle particle distribution function /i (r, p;i), 



1 x/^F 



K 2 3 / 2 



m 



1 + — 



-l 



To (101) 



Ml " their evolution equations follow from the evo- 



,,, ,i , -i , i ji r r lution equation for fi [cf. Eq. (1111) 1. It is ob- 

tells us that n goes to zero and then [cf . M J L ^ V! — V1 

Eqs. m and «] the kinetic coefficients tained the set of evolution equations given in 

a Lo and a n approach zero, i. e. J L (r,t) does Eq. I = for the density, / = 1 for the 

not contribute. first ( vectorial ) flux ' 1 > 2 for the hi S her or " 

der tensorial fluxes, all of which are coupled 

together. 



VII. CONCLUDING REMARKS 

Briefly summarizing the results, it has 
been shown how a statistical nonequilibrium 
ensemble formalism (applicable to the study 
of systems even in conditions far-from equi- 
librium) provides a microscopic foundation 
for a Nonlinear Higher-Order Hydrodynam- 
ics. Its description is based on the set of hy- 
drodynamic variables consisting of the densi- 
ties of energy and matter (particles) and their 

fluxes of all order. In that way we do have a quite gener- 

All these hydrodynamic variables are the alized hydrodynamics under any arbitrary 
average value over the nonequilibrium en- condition of excitation, which, as noticed, 
semble of the corresponding microscopic me- can be referred to as Mesoscopic Hydro- 
chanical operators. Once the complete set Thermodynamics. 
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These generalized hydrodynamic equa- 
tions present on the left side the conserving 
part of the corresponding quantity, and on 
the right-hand side are present the collision 
integrals which include the action of external 
sources and the contributions of scattering 
processes responsible for dissipative effects. 
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Appendix A: Tensorial Coefficients of 

fhk{Q) y ii)2 ( y 2 (2/'-1jii \ O- 

Eqs. dSZD-dSgj and Last Term in Eq. fl62|) ' V ; " V 

The last term on the right of Eq. f l62l) is 

The kinetic tensorial coefficients of Eq. 



(1571) to Eq. (EHD are: 



given by, 



= ^ rfc (Q) (Al) i$(r,f) = + (A9) 



Q 



oo I 



[2fc+2] 



A [2fe+2] = f Tk (Q) Q[ 2fc+2 1, (A2) 

Q 0V/| fc+ '- 1] (r,t)] , (A10) 



oo I 



a? +21 = e /" (o) q m . («> j!i(r ' " = S S p(1 - s) [Ap+21 42 t+ " (r ' ' 



Q 



AM^-^i^r^^fe^) +E(W [Air 2] 0/f + '- 2l (r,t)]}, (All) 
V y fc = 2 

/_-,\fc /A-r«\ fc where the operators p(l, s) and pi are defined 

x 1 ' { 1 1 n[ 2fc+1 l (A4^ 

(2A: — 1)!! V / ' in the main text. 
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Appendix B: The Kinetic Coefficients 
in Eqs. flgl]) and (l62l) 



+ R n (r,t), (C2) 



We do have that 
V 4vr 



T ° (2tt) 3 3 
with k = in Eq. ([II]). And 



y dQ Q A f TO (Q) (Bl) 



5l 2 l (r,t) = -6 Tl /i(r,t)l[ 2 ] + 6a rl -/f(r,t) 



V v^m 2 Q VM 

(B2) 

where ^ (Q) is the Fourier transform of the 
potential energy w (\rj — R M |), n fi is the den- 
sity of particles in the thermal bath, V is the 
volume, and (3~ l = k,BT . Moreover, 



+a L1 - {VI fc (r,t) + [Vl h (r,t)) tr }+R®(r,t), 

(C3) 

where upper index tr stands for transpose, 



M(3 

( 

10 



S®(r,t) = b T A{p 3 [lMl fc (r,f)]} 



777. ^—^ L 



2 1 
— a T 



M/3 7TK 

/dg Q 2 |^(Q)| 2 



(B3) 
(B4) 
(B5) 



5=1 



+ 9a Tl -4 3l (r,t) + 4 3 ](r,t). (C4) 

771 



k j dQ Q 3 \ifj (Q)\ (l + §) 

2 / )7j \ -i 

Kn = -^-^a TO [l + — ) , (B6) 



M/ 



-i 



In which, 

oo 

R„,(r,t)=E{Ai 2fc+2] ©/f +1] (r,t) 

fc=2 

+Af +2| 0V/f](r,t)}, (C5) 



1 /2M/3 
a Li = — \ a TO , (B8J 

OK V 7T 



Appendix C: The Last Terms of Eqs. 
dUD, (fT2|) and d73D 

The contributions present in Eqs. (j7ip . 
and (1731) in section VI are 



R?(r,t)=Z{^o k+2] Ql%%,t) 



k=2 



+1^ (r, t) Af +2] + Af +2] V/f +1 1 (r, f ) 
+ [Af +2 ]0V/f +1 ](r,t)l <r }, (C6) 



T (r, t) = - VV ext (r; *) - Fjvl (r; t) , (CI) 

2 2 
S n (r,t) = 3a ri — l h (r,t)-\ — a L iV/i(r,t) 
m m 



^ 3l (r,t) = ji 3 kr,t) + J^(r,t), (C7) 



[3] 
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+EEp(M)[Af +2, 0r 3, (r,t)]. 

k=2 s=l 

(C9) 

(C8) 



The several kinetic tensorial coefficients 

oo 

jfl (r, f ) = p 3 [A^ 2 ' 1^ (r, t)] , AM are given in Appendix A. 

k=2 



oo 3 



fc=2 s=l 
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